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ABSTRACT 

In this article, we examine the possibility that there exist special scalar-tensor theories 
of gravity with completely nonsingular FRW solutions. Our investigation in fact 
shows that while most probes living in such a Universe never see the singularity, 
gravity waves always do. This is because they couple to both the metric and the 
scalar field, in a way which effectively forces them to move along null geodesies of 
the Einstein conformal frame. Since the metric of the Einstein conformal frame is 
always singular for configurations where matter satisfies the energy conditions, the 
gravity wave world lines are past inextendable beyond the Einstein frame singularity, 
and hence the geometry is still incomplete, and thus singular. We conclude that the 
singularity cannot be entirely removed, but only be made invisible to most, but not 
all, probes in the theory. 
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1 Introduction 



One of the longest-standing challenges to our understanding of gravity has been the 
singularity problem. In any generic theory of gravity and matter, under reasonable 
assumptions about the interactions between particles and fields and about the ways 
of communicating these interactions (which translate into, at least classically, fairly 
loose energy conditions), produces solutions which contain maelstormian regions of 
unbounded curvature In this sense, such solutions seem to point to an intrinsic 
deficiency of the theory that gave them birth, because the very theory that predicted 
such maelstorms looses meaning in these limits. Yet, we must note that not all is 
ill with the fact that such strongly coupled regimes are generic in classical gravita- 
tional theories. Because singular regions generally involve very strong forces between 
particles and very high energies, close to a singularity much of the observed matter 
structure in the Universe can be created starting from arbitrary initial conditions. 
Further, this mechanism is built naturally into the theory, such that extrapolating 
present conditions backwards inevitably results in circumstances under which the 
present could be shaped in a generic way. Indeed, the vast body of astrophysical 
observation does indicate that such a dramatic event, the Big Bang, did take place 
in the past. The cosmological singularities in the past, then, seem perfectly suited 
to encode such cosmogonic furnaces into the theory^. Our task, hence, seems to be 
to bring into accord the facts and the fiction by constructing a theory which could 
encompass all the features of the Big Bang while hopefully not incorporating its own 
demise, in the form of an uncontested singularity. 

While it is not known how, and even if, all singularities may be regulated in a 
theory of gravity, we may be able to gather interesting information by studying the 
difficulties which emerge in the attempt to smooth the edges of the Universe in the 
existing models. A start for any such investigation, of course, should be Einstein's 
theory of General Relativity, since it is in an amazingly good agreement with exper- 
imental observations. In this theory, however, the general theorems due to Hawk- 
ing and Penrose show that the configurations which are determined by the coupled 
Einstein-matter equations of motion, under the assumption that the above-mentioned 
energy conditions are valid, always contain singular regions [|l|. The singularities in 
this context are signaled by geodesic incompleteness. This means that particle trajec- 
tories in such geometries cannot be continued past some hypersurface, because they 
get so strongly focused by gravity's pull that they begin to intersect and hence are 

1 Aside from the cosmological arguments favoring the presence of singularities in any theory of 
gravity, which we describe in the text, there also exist more theoretically-minded arguments, which 
suggest that singularities may be needed in any theory of gravity for consistency reasons. Specifically, 
G. Horowitz and R. Myers have recently proposed that naked singularities in a theory of gravity 
are needed to separate positive mass black hole solutions from the negative mass ones. If the negative 
mass solutions were not singular, they would be more energetically favorable than the empty space, 
thus yielding vacuum instability in any quantum theory one might attempt to build on the original 
classical theory. 
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not smooth any more. This phenomenon does not always imply that the curvature 
of the metric grows beyond bound in such limits. However, the geodesic curvatures 
do, since the geodesic bundles get squeezed very tightly. Hence an observer seeing 
such geodesic curvatures would indeed see a very strong force. With this in mind, the 
notion of geodesic incompleteness is indeed a good indicator of the singularity in the 
space-time. 

The presence of singularities in General Relativity can be taken as a sign that at 
very high energies, or equivalently at very short distances, the theory fails to be a 
completely consistent description of Nature. Because at such high energies quantum 
effects become significant, this suggests that General Relativity should be superseeded 
by the quantum theory of gravity, within which we should seek the answer to the 
conundrum of singularities. Whereas we are still lacking a complete formulation of 
the quantum theory of gravity, there at present is at least one very strong contender, 
string theory. Although the recent developments in string theory, assisted by the 
discovery of the power of duality ||, have greatly improved our understanding of it, 
the theory is still not known in a way that would enable us to ask the questions about 
space-time in a general manner. Instead, we have to either resort to the effective 
action approach which takes into account stringy phenomena in perturbation theory 
, or we could study some special classes of string solutions which can be formulated 
in the nonperturbative regime j5|. The latter approach is clearly more powerful in 
that it allows us to investigate more thoroughly the quantum dynamics of the system 
under scrutiny. But this is available only for some special solutions, most notably the 
BPS states in the string spectrum, and not for any solution we might be interested 
in. In particular, there still does not exist a nonperturbative formulation of generic 
cosmological solutions in string theory. Hence all the investigations of "realistic" string 
cosmologies have been carried out essentially in the effective action approach, which 
is valid for the weak to medium range of couplings and curvatures (see, for example, 
|| - (13|]). We will not dwell on the details of these investigations here, other than to 
mention that the departure of string-theoretic solutions away from General Relativity 
is induced by the presence of additional degrees of freedom which arise in the massless 
string spectrum. These fields, the scalar dilaton field, the torsion tensor (or Kalb- 
Ramond) field, and others, couple to each other and to gravity nonminimally, and can 
influence the dynamics significantly. Thus it makes sense to ask if the spectrum of 
the theory may be tailored in such a way as to produce solutions which do not feature 



any singular behavior - [16|], even if that means abandoning string theory and just 



constructing certain toy models for the purpose of studying the strong coupling limit. 



This approach has been taken recently in (17|, |IBfl . There a special class of scalar- 
tensor theories of gravity was considered, where the matter does not couple to the 
scalar field in some generalized Jordan-Brans-Dicke (JBD) conformal frame, but the 
JBD coupling of the scalar to the metric depends on the value of the scalar field. It 
turns out that by choosing the JBD coupling function one could construct the metric 
in the JBD frame which is smooth over an infinite interval of the JBD comoving 
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time. Since the matter fields couple only to the JBD frame metric, they don't feel 
any singularity at all. However, in the Einstein (E) conformal frame, where the 
metric degrees of freedom have a canonical kinetic term, if the matter sources satisfy 
the energy conditions and consequently the dynamics are still subject to Hawking- 
Penrose theorems, the curvature singularities still plague all the solutions. The E 
frame singularities are invisible to matter probes which do not move along E frame 
geodesies. Rather, in the E frame matter couples to the scalar field, and this coupling 
is adjusted precisely to even out any bumps in the metric. If there were any type of 
probe which couples both to the JBD scalar and the JBD metric, it would detect the 
original singularity. The probe would reach the singularity in a finite extent of its 
world-line, after which it could not be extended any further. This would signal that 
the geometry is still incomplete, at least from the point of view of certain observers. 
Once the subleading interactions are properly taken into account, these observers 
can communicate the presence of the singularity to the remainder. Hence, to check 
whether a solution of some nonminimal theory of gravity is singular or not, it is not 
sufficient to show that curvature invariants are finite in some chosen conformal frame 
and that the space-time geometry is complete in it. Instead, we must investigate 
world-lines of all the physical probes in the theory. Only if all these world-lines are 
not inextendable can the space-time be complete and hence nonsingular. As long as 
there is even a single type of excitation which sees the singularity in any frame, the 
singularity is not absent, but lurks in the geometry waiting to exert its influence on 
the theory. 

Our purpose here is to show that there is at least one such degree of freedom in all 



the scalar-tensor models studied in |17|, [18| . It is the graviton itself. Being generic and 



model-independent, and not immune to the singularity, it resuscitates the singularity 
back into existence of any realistic observer in such theories. To demonstrate this, 
we only need to look at the classical theory. We will present the equations of motion 
of the genuine gravity waves, i.e. the tensor perturbations of the metric, and take 
the geometrical optics limit to show that the wave packets of gravity waves move 
along the E frame null geodesies. These trajectories are past inextendable because of 
the singularity, and the graviton wave packets reach the singular hypersurface after a 
finite extent of the affine parameter along their world-lines. Given this, the singularity 
cannot be completely removed from the geometry: the gravity waves can communicate 
its presence to all other degrees of freedom in a finite time, ultimately making the 
singularity observable. The paper is organized as follows. In the next section, we 
will review the JBD models studied in [O, [18[ , and establish the relationship between 



the JBD and E frames. Section 3. is devoted to the derivation of the gravity wave 
equations of motion, both in the JBD and E frames. The singularity is the central 
notion of section 4. There we will carry out the geometrical optics approximation and 
derive the gravity wave world-lines. In the last section, we will give our conclusions. 
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2 Hiding the Singularity 



Here we will review the models and solutions studied in JT71, |TB| in order to set the 



s = J dVs{x#-— (Vx) 2 - c m (y,vy, g ^)} (2.1) 



stage for out investigation. The theories investigated in these articles were defined by 
the JBD actions with a variable parameter uj{x), given as 

X 

where x is the scalar JBD field, R is the scalar curvature of the JBD frame metric g^ u 
and y and V3^ are any other matter degrees of freedom and their derivatives (field 
strengths). Our signature conventions are g^ v = (— , +, +, +) and R^ u \ a = d\T^ a — .... 
Note that in order to ensure that the gravitational degrees of freedom are not ghost- 
like (i.e. that the graviton propagator never has negative residue), we must require 
X > 0. For simplicity, we choose units such that in the E frame we have 16^^^ = 1. 



The approach taken in [J17 , |18| was to specify the function uj(x) in attempting to 
remove curvature singularities in spatially flat FRW solutions. Note that in the JBD 
frame, the matter fields do not couple to the scalar x a ^ the tree level. To get the 
equations of motion, we can simply vary the action with respect to the independent 



degrees of freedom g^ u , x and 3^- After a bit of straightforward algebra, we find ||T9 



1 _ _ I _o oj /_ I ,_ 2 \ I 



2V,(-V^)+^-^(-)(Vx) 2 = (2.2) 

X ^X X 

V /jT^ Tp,,v 2 ^ ^ u g^iu^m 

where G^ u = R^ u — \g^uR is the Einstein tensor, and is the stress-energy tensor 
of the matter fields y. The conservation of the stress-energy V M T MZ/ = is equivalent 
to the matter equations of motion V At (5£ m /5V At y) = 8£ m /5y. Let us now write the 
explicit form of ( [2.2| ) for the spatially flat FRW cosmologies, assuming that the matter 
stress-energy tensor can be put in the perfect fluid form, T^ u = pg^ u + (p + p)u^u u , 
where is the velocity of the comoving observer, u^u 11 = —1. The FRW metric is 

ds 2 = -dr 2 + a 2 (r)dx 2 (2.3) 

and the comoving velocity is = diag(l,0). The equations of motion (|2.2| ) become 


72 , J 



3H 2 = ^- 3tf * + £ p' + 3H(p + p) = 

2x X 2x 



2H' + 3H 2 + — + 2H— + + — = (2.4) 
X X 2 x 2 2x 

^( x » + 3/fx')+^(- 

X K ' dx K X' 



2 U -{x" + BHjd) + U~)x 2 = 6H' + 12H 2 
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The primes denote derivatives with respect to the JBD frame comoving time r, and 
the JBD frame Hubble parameter is H = a' /a. For simplicity's sake, we will assume 
that the equation of state is p = 7p, with 7 a constant. This is not true in the real 
world, since we know that 7 must be a function of the temperature of the Universe, and 
hence of time. However, since we want to investigate the possibility of constructing 
a classical nonsingular Universe solution, this is a sufficiently good approximation to 
start with. We will determine the limits on 7 later. 



Instead of deriving the equations of motion ( |2.2|) in the JBD frame, we could have 
equally well used the E conformal frame. Conformal transformations of the degrees 
of freedom in a specified theory are merely field redefinitions, and so they cannot 
change physics [M. In this sense, they are to be regarded as changes of reference 
frames, which leave physical observables invariant. The important issue here then is 
clearly to identify the observables correctly. To do that, we must specify an observer, 
building it out of the physical fields in the theory. Once given, this observer doesn't 
care which reference frame we use to compute observables with. So, if we define the 
E frame metric and scalar field by 

9»v = XflV 9 = J d x ^ (2-5) 

(using the overbar to distinguish between the two conformal frames) and, following 
||18|1 , assume that the function t*j(x) is monotonicf] such that the functional relationship 
cf) = 4>(x) is everywhere invertible, we find that in the E frame the effective action is 

S = J d^g{R - (V0) 2 - ^fUy, V3>, ^9r)} ( 2 - 6 ) 

In this conformal frame, the metric kinetic term is canonical, i.e. just R, while now 
the matter fields y couple both to and (p. Note that in order for (|2.5|) to be 
well-defined, we must require 2uj + 3 > 0. This is because when 2uj + 3 = 0, the 
scalar field is not dynamical, since a conformal transformation to the E frame re- 
moves its kinetic term. Further, for 2uj + 3 < 0, the scalar would be ghost-like, 
since its kinetic term would be negative. Now, to find the equations of motion in 
this frame, we can either take the equations (|2.2|) and transform them to the E 
frame using ( |2.5| ), or we can vary the action (|2.6| ). Since in the E frame the mat- 
ter fields couple both to the metric and the scalar 0, the matter equations of mo- 
tion are a little bit more complicated than in ( |2.2| ). Varying ( |2.6| ) with respect to 
y, we obtain V ^[(l/x 2 )5C m /5V ^y] = (l/x 2 )8£ m /5y. Furthermore, recalling that 
5[{l/ X 2 )C m ]/5(j> = -{(dx/#)[2£m + g,u(5C m /5g, u )]}/x 3 = {d X /d<f>)f/2 X , where f 

2 We require that u) is monotonic. But since d(j)/dx — ^f2oj + 3/\/2x> an d we wish that all 
curvature invariants, not only R and W 111 R^ u , are smooth, we should also require that u> is analytic 
|fl8| , because higher derivative invariants depend on the derivatives of 10. We will discuss this in more 
detail following the equation ( [2.1l| ). 
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is the trace of T^,, the resulting equations of motion in the E frame are 

V.T"" = -^TV> = - _ i_^ £m (2.7) 

2x d(j) x og 1 " X 



Again, we will need (|2.7| ) restricted on the spatially flat FRW geometries with perfect 
fluid sources, T^ u = pg^ u + (p + p)u fl u u , where is now the velocity of the E frame 
comoving observer. The E frame metric is 

ds 2 = -dt 2 + a 2 (t)dx 2 (2.8) 

and so u^u^ = —1, u M = diag(l,0). We will show that this is consistent with the 
JBD frame comoving velocity shortly, by establishing the transformation properties 
relating the JBD and E quantities. The equations of motion ( |2.7D then become [TJ| [| 

i? + f + ^ = (2.9) 

+ 3ff0+J-^(3p-p) = O 

Here the dot denotes derivatives with respect to the E frame time t, and the E frame 
Hubble parameter is H = a/a. Using the variational definition of the matter stress- 
energy tensor, we can immediately see that the conformal transformation fl2.5| ) induces 
the change of the stress-energy tensor according to [19] T^ v = T >J, u /x 2 - With this and 



( j2.5|) , we can also easily show that the JBD frame equations of motion (|2.2| ) and 
the E frame equations (|2.7Q map into each other. For the variables describing the 
FRW universes fl2.3|) and (|2.8|) , the transformations are dr = dtf-^/x, cl(t) = a(i)/«/x, 
P = X 2 P an d P = X 2 P- To find out how the comoving velocity transforms, we should 
look at the comoving velocity vector fields in each frame. In the JBD frame, we have 
U — d T and in the E frame U — df, so U — yfxU ■ But this means that the components 
of the comoving velocities in the two frames are identical, m m = = diag(l,0), as 
has been claimed above. This completes our survey of the conformal transformation 
rules for the quantities of interest here. 

Now we can investigate the properties of the theory, using either of the sets of 
field equations ( p.2|) -( p?7[) . We should first recall the Hawking- Penrose singularity 
theorems [T|. To do this, we will use the E frame equations of motion fl2.7p, since the 
gravitational equations of motion are the same as in General Relativity. Consider any 
timelike geodesic with a unit tangent vector field £ M (g^u^^ — — 1) i n a globally- 
hyperbolic space-time (i.e. space-time without any acausal pathologies such as closed 
time-like curves, which we will assume here). The singularity theorems then posit 
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that these geodesies are past-inextendable (i.e. incomplete) as long as the projection 
of the Ricci tensor on the tangent vector field is positive semi-definite, R fJjU ^ fl ^ u > 0. 
By the equations of motion, this condition can be recast as Q l w^ t C + > 0, where 
the tensor M „ is the total stress-energy, M „ = V ^c/>V ' v (f) — (l/2)g^(Vc/>) 2 -\-T tw . This 
requirement is called the strong energy condition (SEC), and is thought to be satisfied 
by most reasonable classical matter sources. In terms of the principal values of the 
stress-energy tensor 0^ = diag(p, ]3i, j3 2 , p 3 ), SEC translates into p + £f =1 ??j > 0, 
P + Pi > 0, i — 1, 2, 3. (Here we use the hat to distinguish the principal values of 
Q^ v from those of T^,, which are both defined in the E frame.) The only feature of 
the JBD system we will be investigating here is the effect of the variation of u with 
X on the singularities, which still needs to be specified. Therefore, we can assume 
that the matter fields y obey SEC, meaning that T lll/ £ t tJ, £ s l/ + T > for all timelike 
geodesies £. For the homogeneous and isotropic cases which we are interested in, this 
tells us that 37 + 1 > 0, or 7 > —1/3. When c/> is included, its stress-energy trivially 
satisfies SEC in the E frame. Hence, also does, being a linear combination of these 
two contributions. We see that in the E frame all timelike geodesies must be past 
inextendable, and thus all cosmological solutions are singular. 

Unfortunately, this does not specify the character of the singularity. Generically, 
geodesic incompleteness as an indicator of the presence of a singularity signals that the 
metric becomes degenerate as some region of space-time is approached. To learn more 
about the actual properties of the singularity, we must look at concrete solutions. If 
we limit our attention to spatially flat FRW cosmologies, which have matter sources 
that satisfy SEC, we will typically find that the singularity arises because at some 
time, say t = 0, the scale factor of the Universe vanishes (or diverges, such as in 
pole-law inflationary solutions) as some power of t, a(t) oc f , « ^ 0[]. Since the 
curvature scalar can be expressed as R = 6H + 12H 2 , we find that in the limit t — > 0, 
R = 6a {2a — 1) / 't 2 + subleading termsfj. So for a 7^ 1/2, the scalar curvature diverges 
at t — 0. If a = 1/2, then near the singularity the Universe is radiation-dominated. 
The scalar curvature vanishes because the radiation sources are conformally invariant 
and hence have vanishing trace of the stress-energy tensor. However, the square of 
the Ricci tensor then diverges, as R^ U R^ U oc 1/t 4 . Moreover, if we look at causal 
geodesies, we can see that in order to move between times t and t = 0, they require 

3 It is clear that hypersurfaces t — const where a(t) is smooth and nonvanishing cannot be singular, 
since the metric there is nondegenerate. We are assuming that singularities can arise only as zeros or 
isolated singularities of a in the functional sense, and further that if a is unbounded for some values 
of t, that such singularities are not essential singularities - i.e. that a is analytic everywhere near the 
singularities, and that if it diverges there it admits a Laurent series expansion with a finite number of 
divergent terms. This is consistent, because essential singularities of a are never a part of the manifold 
in the sense discussed in the text. For assume t = were an essential singularity. Then the geodesic 
distance from the singularity to anywhere else in the manifold satisfies limi_,o AA oc Ym\ t ^Q ta(i) . 
But since t = is an essential singularity, lim t ^o ta — > 00. Hence, AA always diverges. 

4 Note that this is valid even for the cases when the scale factor vanishes faster than the power 
law, such as a oc t 2 / 3 (ln(t)) 1 / 2 , which is known to be the limiting form of the solution dominated by 
composite hadrons pl| . 
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a lapse of the affine parameter equal to 



AA= / dt a/Vv 2 + m 2 a 2 (2.10) 



o 



which can be obtained by solving the causal geodesic equations for the metric (|2.8|) , 
and where v 2 and m 2 arise as constants of integration; v 2 is a nonnegative constant, 
and m 2 = 0, 1. In the limit of small t, when a oc t a , we can approximate this expression 
with AA oc t 1+a / y/v 2 + m 2 t 2a , hence noting that we can always choose v 2 and m 2 such 
that AA is finite^]. This means that there always exist causal geodesies which reach 
out of the singularity to anywhere in the Universe in a finite proper time - i.e. they 
are past-inextendable and so incomplete. As a result, all the FRW solutions with this 
type of behavior are indeed singular. 

To see how the mechanism of conformal transformations could regulate the sin- 
gularities, recall that since all FRW solutions (that is, k = ±1,0) are conformal to 
static geometries with maximally symmetric spatial slices, the Weyl tensor of these 
solutions is a constant and hence does not encode any information about the Big Bang 
singularity. The information about the singularity is completely encoded in the Ricci 
sector of the curvature, which changes under conformal transformations. Indeed, we 
can look at the leading form of the curvature in the vicinity of the singularity of 
every solution we have discussed above. Applying the field redefinitions (|2.5|) to the 
solutions, and using the equations of motion in the E frame fl2.7|) , we can show that 
the JBD and E frame curvatures are related by 

*^ +i ^-Sf(^ + 3)} (,n, 

We have seen above that as the singularity is approached, the E frame Ricci scalar 
diverges as R oc 1/t 2 . Hence if in this limit x ~ * t with (3 > 2, the divergent 
contribution of R in (|2.11 ) can be tamed. What about the other two terms? It is 



straightforward to see that in the limit a(t) — > t a , the equations of motion (|2.9| ) are 
approximated by 

a 2 _ 2 p a _ 2 1 + 7- 

t 2 ~ 2 + 2 t 2 ~ 2 + 4 9 

OL OL • A 

p + 3(l + 7 )-p-A<j>p = + 3 -0+-p = O (2.12) 



where A = (3 / y—l)/^Auj(xo) + 6 > when 7 > —1/3 is a finite constant since 2lu+3 > 
0. On the other hand, we can see from Q2.2|) that we can always write down the exact 
solution for the fluid, since in the JBD frame the fluid couples only to the metric. The 
solution is p = po/ a - 3< " 1+7 ' 1 - Now, since p = p/x 2 , we can transform the energy density 

5 e.g. for a > — 1, take m 2 = and v 2 = 1; then, AA oc t 1 ^ 
(static observers!), so that AA oc t. 
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to the E frame: p = p X ( ' 37 ~ 1 ^ 2 /^ 3< ' 1+7 ^ m the limit t — > 0, as a oc t a and x ^ 
we find p oc £( 3 7-i)/3/2-3a(i+7) Hence in order for the p pole contribution to the JBD 
curvature to be smoothed, we must require /3+ (37 — l)/3/2 — 3a (1 + 7) > 0. From the 
equations of motion near the singularity, we see that p < 6a 2 /t 2 for all t. Therefore, 
for all solutions, (37 — l)/3/2 — 3a(l + 7) > —2, which means that the regulator 
field x is 110 worse than in the previous case. Let us now compare this to the scalar 
field contributions. Assuming a strictly greater than order relation in the previous 
inequality, so that p falls off slower than 1/t 2 , we see from the constraint equation 
solved for 0, 2 /2 = 3H 2 - p, that near the singularity the dominant contribution 
comes from H. This means that the scalar field is oc y/Ea/t. But then, sufficiently 
close to t = 0, the p-dependent contribution to the equation is also negligible. 
As a result, in this limit a — > 1/3, and hence the solution approaches a — > t 1 / 3 , 



+ y2/31n(t) with negligible matter field contributions - i.e. we get the 
scalar field-dominated cosmology [§, [Kj. Inspecting (|2.12 ) near the singularity we 



can see that the matter sources can never dominate over the scalar field[]. Hence near 
the singularity, the 0-dependent contributions can never be subleading to the matter 
stress-energy contribution, regardless of the kind of matter field. 

The only remaining possibility is that the two sources remain of equal importance 
near the singularity. Indeed, a glance at ( |2.12j ) suggests that a solution of the asymp- 



totic form oc p oc 1/t 2 near the singularity is admissible. Yet, a closer look 
reveals that, unless 7 = — 1, the four equations in ( p. 12 ) are consistent only if p = 



- hence again leading to the scalar field-dominated solution. If 7 = — 1 (which corre- 
sponds to the JBD frame cosmo logical constant, since p = const from the equations 
of motion (|2.4j) ), the solution which treats the sources in an egalitarian way is admis- 
sible. However, in this case the parameter (3 cannot be adjusted to be > 2. It is fixed 
by 7 = — 1 and 3a (1 + 7) — /5(37 — 1)/2 = 2 to be f3 — 1. So the JBD frame solution is 
still singular! For this case, it is impossible to remove the E frame singularity by going 
to the JBD frame via the conformal transformation ( j2.5| ), regardless of the form of 
the coupling function w(x)0. From this discussion, we can see that all the JBD frame 
solutions without a singularity must have a unique behavior near the singularity in the 
E frame - that of the scalar field-dominated cosmology. In fact, this is the reason why 
a conformal transformation can be used to smooth the solution. Once the solution is 



6 Suppose they did; then, near the singularity the scale factor would behave as a a £ 2 /3(i+7) 
and hence p oc 1/t 2 . But this then produces the response in the scalar field according to <j> = 
P/t 2 / 1+1 + Q/t, where P and Q 7^ are integration constants. 

7 For completeness, we should also mention another special example where frame switching fails, 
although this case does not violate the results of |l7], since the scalar field is constant. Namely, 
if the matter is in the form of pure radiation, there exist solutions with <p = const, as can be seen 



from either (2.2) or ( |2.7| ). Hence the JBD and the E frame are identical, and so the singularity is 
not removed. Nevertheless, these solutions are not past attractors for all late-time radiation-filled 
Universes, since as long as <f> ^ 0, early enough it will dominate over radiation. So these solutions 
are not generic, but they illustrate that scalar field must dominate near the singularity in order for 
the frame switching to be successful in smoothing the solution. 



9 



dominated by the scalar field near the singularity, we can choose an almost arbitrary 
function of the scalar field to conformally transform the solution with - as long as the 
conformal factor vanishes in the limit t — > at least as fast as t 2 . 

To complete the smoothing of the singularity within the context of the JBD theory, 
we must require that the hypersurface t = can never be reached by any particle in 
the JBD world. If we look at the matter sector of the theory, in which fields move 
along the JBD frame geodesies, we should require that no such geodesic ever reaches 
the t = hypersurface. When we solve the geodesic equations of the JBD metric 
( |2.3| ), and express them in terms of the E frame time for convenience, we find that 
the JBD frame geodesic lapse between hypersurfaces ti and t is 



where as before v 2 > and m 2 = 0, 1 are integration constants. To make sure that 
all the geodesies are complete, we must require that AA diverges as either of the 
limits of integration goes to zero. Substituting the asymptotic form of the functions 
a and \ which we have deduced above, we find that as long as (3 > max(2, a + 1), 
the geodesic distance AA between the t = hypersurface and any other t = const 
hypersurface diverges. Further, if we set v 2 = 0, m 2 = 1, we recover the integrated 
coordinate transformation between the JBD comoving time r and the E comoving 
time t, t = J l dt/x, as we should. So we see that the hypersurface t = as seen 
from the JBD frame corresponds to the infinite past (or future!) of the solution, as 
long as (3 > 2: lim^o r ~~ > ±oo. Therefore, the conformal transformation to the JBD 
frame does push away the singularity to an infinite distance, while making the JBD 
curvature finite. 

At this point, we must remember that while we are using the conformal transfor- 
mation (|2.5| ) to regulate the solution at t — 0, we must make sure that it does not 
behave badly elsewhere. The conditions which oo(4>) must satisfy in order to guarantee 
that the JBD solutions are nonsingular can be summarized in terms of Q = 2uj + 3, 
following [)I8| , as 

(1) d n Q/dx n are smooth everywhere except in the limit \ — > 1 for all n G Z + , 

(2) as x "~ 1) both Q and it's derivatives diverge as some inverse power of (1 — \)i 
and 

(3) finally and most importantly, that near the singularity the solutions are dominated 
by the JBD scalar field <fi (or equivalently, x), which implies lim x ^o^ > 1/3- 

We emphasize here that whereas this condition is vacuous for a range of admissible 
types of matter, it is not automatically true for all solutions. In particular, we have 
given two simple examples (with the JBD frame matter being either the cosmological 
constant or a radiation fluid) where the scalar field in the E frame is not dominant 
near the singularity and hence the JBD frame solutions remain singular. 




(2.13) 
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However, as we have indicated earlier in the discussion preceding the investigation 
of the conformal pictures of the physics given by (|2.1|) - (|2.6|) , even if the conformal 
rescaling did produce a smooth metric in the JBD frame, it may not have really re- 
moved the singularity if there still remained even a single degree of freedom which 
coupled to the E frame geometry. In what follows we are going to show that gravita- 
tional waves move along worldlines in the JBD frame which are not geodesies. Rather, 
they are deformed by the dilaton force such that they are identical to the E frame 
geodesies! This should not be a complete surprise since it is the E frame where the 
graviton kinetic term takes the canonical form, being just the Ricci scalar. It has 
been argued by Shapere, Trivedi, and Wilczek that in theories of antisymmet- 
ric p-forms coupled to scalars via f(4>)F^ l 4l+l the tensor field quanta move along 
geodesies of the metric where the term f((p) is absorbed away by a conformal trans- 
formation. This was later proved in |23[] by the present authors and collaborators 
to hold as the geometrical optics limit of the p-form field equations. Since gravity 
is a gauge theory like the p-form field theories, when we apply the same technique 
to the equations of motion of gravitational perturbations in the FRW background, 
we find that the gravity waves in general see a scalar force. Therefore since the E 
frame metric is singular, the gravity wave worldlines are past-inextendable, and hence 
the full physical arena of the theory given in either frame is still incomplete! This 
must be taken as a signature of a latent singularity, which simply didn't go away by 
a conformal transformation, but only appeared invisible to the matter sector probes. 



3 Gravity's Redoubt 

Below we will consider small perturbations away from a fixed curved background, 
and look at their dynamics to the lowest order in the fluctuation. Since we want to 
study only the pure gravitational excitations, we will assume that the matter and 
the JBD scalar are unperturbed, and impose the transverse traceless gauge on the 
metric fluctuation. Further, since we will look for the perturbations around the FRW 
solutions, we will impose the stationary gauge on the background metric, which will 
therefore leave the background solutions identical to the ones studied in section 2. 
The resulting equations of motion give the correct dynamics of the two independent 
graviton polarizations, the + and the x modes, in exactly the same way as discussed 
previously EH]. In the E frame, these equations will be the same as in the ordinary 
General Relativity, because of the gauge conditions. In the JBD frame, they will 
contain an additional coupling to the scalar x, which arises from the noncanonical 
form of the graviton propagator. 

In the E frame, therefore, we will use the following ansatz for the metric and 
matter fields: 

V = fa + V = W = (3. 14) 
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The genuine gravitational degrees of freedom correspond to two graviton polarizations 
defined by imposing the transverse traceless conditions on the metric perturbation. 
In the momentum space, this means that we require the polarization tensor of the 
graviton e pv to be traceless and orthogonal to the direction of motion, i.e. V 1 ^ = 
k^e^u = 0. In terms of the perturbation h pv , these conditions can be written down as 
h^n = Vfji^v = 0. Here the covariant derivatives and raising and lowering of indices 
are taken with respect to the background metrics g^ iu . The inverse metric, to linear 
order in perturbation h is g^ u = g^ v — . Notice here that the determinant of the 
metric is not perturbed: g = det(g „„) = g + p — g. A straightforward calculation 

shows that the Christoffel symbols in the perturbed background are T uX = T^ x + 7^ A 
where 7^ A = (l/2)g w (V \h up + V u h pX — V p h v \). The Ricci tensor, expanded to linear 
order in h pu , is, in the transverse traceless gauge, 



Tl^u = R^u + ^V A V M /i\ + ^V x V u h\ - ^V 2 V (3.15) 

Note that the perturbation in R pv is traceless because of the radiation gauge con- 
ditions, as expected. Tracing out the gravitational equations of motion in ( |2.7j ) and 
subtracting the Ricci scalar-dependent term, we see that R pv = V^V^ + T pu — 
(l/2)g fMU T. The perturbed equations of motion are identical to these. In order to find 
the equations of motion of the perturbations, then, we have to expand the equations 
to linear order in and cancel the lowest order terms since the background is a so- 
lution to (|2.7| ). We also need to demonstrate that the perturbation is consistent with 
the matter sector of the theory, in that it is not exciting any perturbations there. 
At this time, it is convenient to refer to the explicit form of the background solu- 
tion, and introduce the gauge conditions for it. We will work with the synchronous 
gauge, where the metric is exactly the same as in ( |2.8| ), ds 2 = —dt 2 + a 2 (t)<ir 2 , and 
the perturbation satisfies h 00 = h ok = 0. Then, on all cosmological backgrounds of 
the previous section, h^V u (p = 4>h^° = 0, and so, the perturbed d'Alembertian is 
identical to the unperturbed one: 

^2 , 1 „ , /r^„ „ _ , 1 



V = -fzdMggrW) = V 2 - —d^h^d^) = V 2 (3.16) 

Next, since we are working in the synchronous gauge, the comoving velocity is un- 
changed: v! 1 = = diag(l,0), and so is the trace of the stress energy tensor. The 
perturbation of the stress-energy tensor for fluid sources is 

vT P v = T p u ~ ~^Q P vT — T pv + —g^JT 

= phfiu - Ifi^T (3.17) 

- p+P 
^ 2 

Also, the stress-energy conservation equation is identical to the case without the per- 
turbation because of the gauge conditions; specifically, V^T 7 "' = V^T^. Substituting 
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these conditions and the equation ( |3.15| ) in { \2.7\ ), we obtain the covariant form of the 
equations of motion for tensor perturbations: 



V A V^\ + V x V u h\ = V% u + V (3.18) 

We now want to rewrite this equation in terms of the mixed perturbation tensor h^v, 
which is the natural variable to use because of the gauge conditions. Using the explicit 
form of the background metric ( |2.8|) , we can verify that in the synchronous gauge the 
following conditions hold identically (see, e.g. 



Vy/ioo — Vo/io — Vj/ifeo — —Hhj k 

Vo/ijfe = hjk — 2Hhjk Vj/ijfc = dihjk (3.19) 

Since these equations are covariant with respect to the assumed background, we can 
raise the indices using the background metric, and after some straightforward algebra, 
we can also verify the following conditions on the second covariant derivatives of the 
perturbation []24], [25| : 

V 2 P = V 2 /^- = V m V°^q = V^V /^- = V M V^ = 

V fl V r W k = (B + 4H 2 )W k (3.20) 

W, = -X 2 h^, - Ji\ - 3H%\ + 2E 2 h\ 
a 1 

In these equations, V 2 is just the three-dimensional flat space Laplacian, V 2 = £:- =1 <9 2 . 



Upon substituting these expressions in ( |3.18| ), we find 



Ji\ + mti\ - \v 2 h j k + (2H + QH 2 + V - 1 /-)h? k = (3.21) 
a z 2 

The last term looks like the environment-induced mass; however, by the E frame FRW 
equations of motion (|2.9|) , this term is identically zero. Hence, finally, the equations 
of motion for the transverse traceless perturbations of the metric are 

fi\ + 3Hh\ - \v 2 hh = (3.22) 

cr 

i.e. the propagation equations for a set of minimally coupled scalars. The index 
structure of the perturbations h? k can be easily accounted for by going to the mode 
expansion h? k = P k (t,p) exp(ip-x). The gauge conditions e 7 j = PjP k = can be easily 
solved as follows. We orient the spatial reference frame s.t. the z axis is along the 
direction of propagation of the wave, and then we find that two linearly independent 
polarization tensors are given by two Pauli matrices e + = 03 and e x =01. Any other 
perturbation is their linear combination: e j k = f+e + + f x £x, where f k , k = (+, x) 
are the mode functions. Since we will see later that the mixed index perturbations 
h^ v are conformally invariant, and f k are modes of these degrees of freedom, f k = f k 
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and from now on we will omit the bars from f^s. In a different coordinate system, 
the basis polarizations are given by = R~ 1 {p)<JkR{p) ) where R(p) is the rotation 
matrix which orients p along the z axis. Then, as we mentioned above, the equations 



of motion (|3.22|) can be rewritten in terms of the mode functions /+, f x as the Klein- 
Gordon equations for a set of minimally coupled massless scalar fields, exactly as in 
General Relativity: 

V 2 /fc = (3.23) 

This equation suggests very strongly that the metric perturbations choose to propa- 
gate along geodesies in the Einstein conformal frame. We will demonstrate that this 
is indeed true in the next section. 

Now we turn our attention to the description of the wave propagation from the 
point of view of the JBD frame. We could proceed in precisely the same way as in the 
derivation of the wave equation in the E frame. Starting with a fixed JBD background, 
we add a perturbation to the metric as g^ u = g^ v + h^ u , and expand the equations 
of motion ( |2.2j ) to first order in the fluctuation. Here, however, we don't have to 
repeat all the steps of the derivation in the E frame, because most of the details are 
the same, with barred quantities (E frame) replaced with the unbarred ones. We will 
therefore outline only the main points and establish a correspondence between the two 
pictures. First, note that since the field redefinitions ( |2.5|) imply that = h^/x, the 
mixed index perturbations are conformally invariant: h fl u = h^ v . Next, the transverse 
traceless conditions in the JBD synchronous gauge are also conformally invariant: 
clearly, the equality of the mixed index tensors implies that if one is traceless, so is 
the other. Further, the constraints of the JBD synchronous gauge on hF v are also 
identical: h°o = h°j = 0. Finally, since 



= V^ v -2h%^+h%^ = VJi^ (3.24) 



the transversality condition is also the same: V^/i^j, = 0. Hence the perturbation 
of the Ricci tensor in the JBD frame is of the same form as (|3.15|) , except that the 
quantities are all unbarred. There however are additional source terms, since now we 
have to take into account the contributions of the second derivatives of the field x as 
presented in ( |2.2| ). These terms will be important for showing that the geometrical 
optics limit of the wave dynamics picks the E frame evolution. Using the conformal 
correspondence we have established so far, we can write down the equations of motion 
for the perturbations in the JBD frame: 



X 

^(VV-V^-V^;) (3.25) 
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In the FRW background, after some straightforward algebra, we can rewrite these as: 



h"\ + 3Hh'\ - ^V 2 U k + ^h'\ + (2H' + 6H 2 + X " + 5Hx> + V —^)h\ = (3.26) 
a X X 2x ' 

Again, the environment-induced mass term vanishes, now by way of the JBD frame 



FRW equations of motion (2^). The final set of the equations of motion for the 



transverse traceless perturbations of the metric are 

h"\ + 3Hh'\ - \v 2 V k + ^h'\ = (3.27) 
a 1 x 

These equations contain the term proportional to x'/Xi since in the JBD frame the 
tensor perturbations also couple to the scalar \. To see how this term arises in the 
form above, recall that because ( 3.22|) and ( |3.27 ) map into each other under confor- 



mal transformations ( |2.5|) , the x'-dependent term must arise from the effect of the 
conformal map on the comoving time. Now, it is evident that the same polarization 
basis as that used in the E frame analysis above can be employed to represent an 
arbitrary perturbation. In fact, we could just take any solution of p. 22] ), and confor- 
mally transform it to the JBD frame according to (|2.5|) , and it will be guaranteed to 
solve ( |3.27| ) as well. In terms of the mode functions, then, the JBD frame equation of 
motion for gravity waves (|3.27|) can be written as 



VJxV/fe) =0 (3.28) 



where f k are the mode functions, with k G {+, x}. In this case, the modes propagate 
under the influence of an additional conformal coupling to \. This equation is confor- 
mal to ( |3.23| ) under (|2.5|). These two equations, ( |3.23| ) and ( |3.28| ), will be our starting 
point of derivation of the geometrical optics approximation in the next section. 



4 Singularity Revealed 

We are finally ready to show that smoothing the JBD frame is not sufficient to remove 
the singularity from FRW cosmological solutions we have discussed in section 2. As 
we have indicated in the introduction, this is because the gravitons, when considered 
as probes of the geometry in the classical limit, move along null geodesies of the E 
frame, and not JBD. Since the E frame is manifestly singular, the graviton worldlines 
are incomplete and so they reintroduce the singularity's effects even in the JBD frame. 

Let us first derive the geodesic equations and then consider their implications. 
As a warm up, we will first find the geodesies in the E frame, since there the mode 
equation is very simple - it is just the Klein-Gordon equation for a minimally coupled 
massless scalar field for each polarization mode of the graviton sector, wfk = 0. Now, 
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the geometrical optics limit corresponds to setting f k = exp(iS), and identifying the 
phase S as the action of the pointlike probe which replaces the wave packet |23| . The 



field equation in terms of S becomes iV 2 S — (VS) 2 = 0, and so, assuming that S is 
real and separating the real and imaginary parts of the equation we find 

(VS) 2 = V 2 S = (4.29) 

Now, after S is integrated, it must be a function of the space-time coordinates only, 
since it is the phase of fk- Hence, if we use S = S(x^), the action must be representable 
as a path integral along a geodesic which the wave packet is following: S — j dx^V^S. 
Let us now introduce = V^S and recall that dx^ 1 = x^dX along a geodesic, where 
A is the afline parameter. Hence [^J, 

3= f V^dX (4.30) 



The equations of motion for S ( |4.29| ) become, in terms of the field V^, the following two 
constraints: g^V^Vv = 0, V^g^V^ = 0. The first constraint is local and hence can be 
easily enforced at the level of the particle action with the help of a Lagrange multiplier. 
It just tells us that is a null vector. The second constraint is not local, and in fact in 
the geometrical optics limit is always a very small quantity compared to the first one, 
and is usually ignored. Here we will retain it, and use it to determine the Lagrange 
multiplier. The constrained particle action is then S = J dXiV^x^ + fjg^V^Vv) . In the 
geometrical optics approximation, because of destructive interference of waves, only 
those trajectories for which S is extremized survive. Treating as an independent 
variable, and varying 5* with respect to it, we find g flu V u = —x^/2fj. This simply means 
that Vfj, is tangent to the wave packets world-lines, which by g^V^Vy = must be 
null. The differential (second) constraint translates into the condition V^x^/fj) = 0. 
Now, in the flat space limit, the covariant derivative would have become an ordinary 
derivative, and we would have been able to use d^x^ = c?(9 /J x M ) /dX = to assert 
that x^d^l/fj) = d{l/fj)/dX = 0, i.e. that fj is a constant along trajectories. Then, 
with an appropriate normalization (fj = —1/4) we would have found S = J dXx^x^ 
with x^x^ = 0, i.e. just the standard action of a relativistic massless particle. In 
the curved space, the tangent vector field of a congruence of geodesies need not be 
divergenceless in general, i.e. it need not be V ti x fl = 0. However, since this quantity 
is the second constraint in ( |4.29| ), because of destructive interference, its contribution 



to wave propagation is negligible in the geometrical optics approximation. Moreover, 
if we go to the Riemann normal coordinates, defined at a point x by the condition 
f^ A (x ) = 0, we will indeed find that to the lowest order V^i^ vanishes. So, we will 
set V^i^ = 0, as a part of the geometrical optics approximation. This then tells us 
that x^V^l/fj) = d(l/fj)/dX = - just as in the flat space limit, fj is a constant along 
geodesies. Therefore, the particle action in the E frame becomes, after again choosing 
fj = -1/4, @ 

S = f dXg^x^x" g^x" = (4.31) 
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which defines dynamics of the E frame massless minimally coupled particles - i.e. the 
null geodesies. The geodesic equations can be derived straightforwardly by varying 
this action and imposing g fJLV x^ 1 x v = 0. The result is 

x^ + f l x x v x x = g^afa? = (4.32) 

Hence as claimed, in the geometrical optics limit the quanta of gravitational pertur- 
bations move along the E frame null geodesies, and not along trajectories in the JBD 
frame. 

We can now derive the same result using the JBD form of the equation of motion, 
in order to check how the dilaton force deforms the JBD graviton trajectories into the 



E frame null geodesies. To show that the result is a conformal transform of ( 4.32 ), we 
will derive it from first principles, rather than just apply the field redefinition ( |2.5| ). 
Bearing in mind that we are following the motion of the same mode {h^ v = /i M v ), but 
in a different reference frame, we again use fk = exp(iS), but now we substitute it in 
the JBD frame equation V^(xV M /fc) = from Eq. ( |3.28|) . When S is real, we again 
get two equations, 

(yS) 2 = V M (xV^) = (4.33) 

The first of these two equations is identical, up to an overall factor of Xi t° the first 
equation in ( f4.29|) . Thus it also must correspond to the null condition, g flu V fl V u = 
where we still have = V^S, albeit in the JBD frame. The second equation is a 



little bit more complicated than in ( |4.29| ). To see the effect of the coupling Xi let us 



rederive the equations of motion for graviton probes along the same lines we followed 
in the E frame calculation above. As we have shown, S is a function of the space-time 
coordinates only, and when we use S = S(x fl ), S again is just the (same) integral 
S = J dx^V^S. This of course must be identical, since V^S* = d^S = V^S. With 
Vfj, = V^S 1 and dx^ 1 = x^dX along a trajectory, 

S = J V^dX (4.34) 

The JBD frame equations of motion for S then are g^V^Vy = 0, V M (x^ /iV V l/ ) = 0. 
Enforcing the first constraint at the level of the action with the help of a Lagrange 
multiplier, we find S = J dX(V^x^ + rjg^ u V^V u ). Note that here we are contracting 
the indices with the JBD frame metric g^, and so rj ^ fj. Varying this action with 
respect to V^, we find g^ v V v = — x^/2^, and so again, is a null vector tangent 
to the graviton trajectories in the JBD frame. However, the differential constraint 
is then V M (xi M /^) = 0. We still implement V^x^ = as a part of the geometrical 
optics approximation, and find x^V ^ix/v) = d(x/v)/dX = 0, implying that x/v = c 
is constant along each particle path. Choosing c = —4, we can write the JBD frame 
particle action as 

S= f dX-g^x u g^x u = (4.35) 

J X 
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The particle trajectories are those paths which extremize this action. By varying this 
action, and taking the null constraint into account, we find p3| 



x M + T^ x x u x x - —x v x» = g^xPi? = (4.36) 

Hence in the JBD frame gravitons do not move along geodesies, but rather along null 
trajectories determined by the additional force proportional to the 4-gradient of the 
scalar field x- Of course, this is what we have expected all along, since we can see 



that in the JBD frame action (|2.1| ) the \R coupling implies that the JBD "gravitons'' 



have x field charge, and hence must couple to x's field strength. Yet, when we use the 



conformal transformation (2.5) which transforms the JBD action (|2.1|) to the E frame 
action ( |2.6|) , rendering the graviton kinetic term canonical, i.e. just R, we see that 
this same redefinition removes completely the x-dependent force from the equation 



( (4.36| ). Under the conformal transformation ( |2.5| ), the connexion changes according 



to = F^ A + (1/2x)(<5^VaX + S^V u x — g^QvX^ pX) an d when contracted with x v x x , 
the difference precisely cancels the x-dependent force in ( [4.36| ). Hence, the equations 
( |4.32|) and ( |4.36| ) are conformal images of each other, and they imply that the gravity 
wave packets, to the lowest order, move along null geodesies of the Einstein frame 
metric g^ v . 

Having proven this, we conclude that the conformal removal of the singularity was 
only partially successful. By the construction of the theory as given in the action 
Q2.1| ), the matter fields y coupled only to the JBD frame metric g^. This metric was 
made smooth by the choice of the coupling function o>(x), whose role was to push the 
E frame singularity to the asymptotic timelike infinity of the geometry. Due to this 
effect, the point-like matter probes were unable to ever reach the singular region. Since 
they follow the JBD frame geodesies, these geodesies are all complete. However, in this 
case, geodesic completeness is insufficient to conclude that the solution is nonsingular. 
The gravity waves couple to both the metric and the scalar field in the JBD frame, and 
the net effect of these couplings is to deform the graviton trajectories back to the E 
frame null geodesies. These worldlines are incomplete, because the E frame metric has 
a singularity a finite affine distance from any other place in the manifold. As a result, 
we cannot arbitrarily extend the history of one such Universe. Sooner or later, we will 
reach the E frame singularity, where we will have to deal with the problem of defining 
the proper initial conditions for gravity waves. Because the graviton sector there is 
ill-defined, we simply wouldn't be able to unambiguously set the initial conditions. 
Gravity waves could then communicate the presence of the singularity to the other 
degrees of freedom in the theory, at higher order. To see this, note that because the 
matter fields couple to the metric, they would also couple to the tensor perturbations 
at higher order. Since to linear order the tensor perturbations were the only degrees 
of freedom which detected the singularity, none of the other modes could be adjusted 
to completely cancel the influence of the singularity on the gravitons. Therefore the 
gravitons would render the hidden singularity again visible to all the fields in the 
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theory. In particular, the relic gravitons from a very early era of the Universe that 
survived recombination and that comprise the present graviton background would still 
encode information about the maelstorm they came from. 

5 Conclusion 

In this article, we have shown that conformal transformations alone cannot completely 
remove the initial singularity from a cosmological solution in the scalar-tensor models, 
even if we allow the scalar-gravity coupling strength to depend on the scalar field. The 
matter sector in the models we have considered consists only of modes whose stress- 
energy tensor satisfies the strong energy condition. Clearly, if the matter sources 
violate the SEC, then evading singularity may yet be possible - but this cannot be 
accomplished by a simple conformal transformation involving the JBD scalar. Such 
conformal transformations merely hide the singularity from observers which propa- 
gate along geodesies of the JBD frame metric, that can be arranged to be smooth 
by adjusting the coupling u(<j>). However, once different observers are allowed, which 
couple to both the JBD frame metric and scalar, because of the additional scalar force 
they don't move along the JBD frame geodesies, but along those of a different con- 
formal frame. In particular, gravitons move along the Einstein frame null geodesies, 
which are always incomplete, because the E frame solutions are all singular. Hence 
gravity waves see the singularity in all scalar-tensor models, regardless of the specifics 
of the model in question. We can in fact see that a similar property should hold 
in any nonminimal effective theory of gravity. The plain vanilla field redefinitions 
cannot remove singularities because they do not couple universally to all the fields in 
the theory, and so there will always be a mode which will discern the presence of the 
singularity in the manifold. 

In string theory, however, this may merely mean that the effective action approach 
must break down close to the singularity. Since all the states in string theory are 
comprised of strings (and D-branes, as we have seen recently in some of the models), 
close to the singularity, gravity's pull on low energy point-like degrees of freedom is so 
strong that they effectively decompactify. Thus near the singularity instead of a gas 
of interacting highly energetic particles we find a gas of interacting highly energetic 
strings. But strings couple naturally to the string frame metric, which therefore 
seems to be the frame we must choose to study the effects of finite size of probes (or 
higher order a' terms in the derivative expansion). However, strings also couple to the 
dilaton, and we must consider its effect on dynamics too, since it represents the string 
coupling constant and so controls the validity of the semiclassical approximation. 
Therefore in string theory, if we are to remove the singularity, we must regulate both 
the metric and the dilaton, which is reminiscent of the situation in the scalar-tensor 
models we have considered here. This also implies that if we manage to cure the 
solution in one frame, as long as we take only reasonable field redefinitions (i.e. those 
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which don't alter the global properties of the space-time) we find that the solution is 
free of singularities in any frame. This, of course, is consistent with the fact that field 
redefinitions do not change physics, but just alter the language we describe physics 
with. 

Finding whether the singularities could be removed from the effective theory still 
awaits. What we have learned to date seems to indicate that to regulate the singulari- 
ties we must impose severe alterations on the theory. Simple tweaks of the classical (or 
semi-classical) effective action do not seem to work. To illustrate this point, we may 
recall the graceful exit problem in the Pre-Big-Bang scenario |TJj]. The idea, briefly, 
was to use the multiplicity of solutions which arise because of duality, and paste them 
together in such a way that the result is a smooth Universe of an infinite life span 
and with a region of very large coupling and curvature, mimicking the Big Bang. The 
difficulty with this proposal is that the continuity of the solutions and the equations of 
motion, in the effective potential approximation (i.e. the effective action truncated to 
second order in derivatives) makes the smooth matching of the two branches with the 
proper asymptotic behavior impossible unless very exotic conditions develop at very 
high energies ||16|| . The matching is supposed to take place very near the singularity 



- but in the semiclassical approximation, the singularity does all the steering of the 
dynamics near it, and hence the solutions cannot evade it. The branch changing thus 
is not possible in the semiclassical limit, with the matter sources satisfying an even 
weaker version of the energy conditions than in General Relativity. A necessary condi- 
tion (but not sufficient) for branch changing has been derived recently by R. Brustein 



and R. Madden (in the last of Ref. [fLq|), who showed that unless the null energy 
condition (i.e. the requirement p + p > for fluid sources) is violated, the branch 
change cannot occur. This still does not guarantee the evasion of the singularity via 
branch-changing, but it does tell us that unless we violate the null energy condition in 
some way, in a region of high curvature, we cannot even hope to avoid the singularity. 
Since the conventional matter (point-like or string-like) degrees of freedom generally 
do not violate the null energy condition, this requires the presence of exotic matter 



sources ("string phase", [p6f] ) or nonperturbative phenomena (such as discussed by 



T2] , |T5T , to account for singularity smoothing. At this moment, we know very little 
about this type of matter. The recent rapid development of string theory however 
gives hope that we may learn more in foreseeable future! 
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